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In the present work, gravitational collapse of an inhomogeneous spherical star model, consisting 
of inhomogeneous dust fluid (dark matter) in the background of dark energy is considered. The 
collapsing process is examined first separately for both dark matter and dark energy and then 
under the combined effect of dark matter and dark energy with or without interaction. The dark 
energy is considered in the form of perfect fluid and both marginally and non-marginally bound 
cases are considered for the collapsing model. Finally dark energy in the form of anisotropic fluid 
is investigated and it is found to be similar to ref. [12]. 
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I. INTRODUCTION 



In recent past, there are remarkable observational evidences which contradict the present day prediction of 
Standard Cosmology. The observed data from the high red shift of type I-a supernova [1] suggests that the 
universe at present is accelerating instead of deceleration (a prediction of Standard Cosmology). This was 
confirmed by the observed data from measurements of the fluctuations in the power spectrum of the cosmic 
, microwave background radiation [2] and large scale structure [3] . 

\ To incorporate such accelerating phase within the frame-work of Einstein's general relativity one requires 
source of repulsive gravity, termed as dark energy. This component of the matter distribution of the universe 

\Q ' should have a large negative pressure and hence violates the strong energy condition and dominates over matter 
at present. Also astronomical observations predict that at present the universe contains approximately | dark 

""^ | energy and | dark matter. Although, the nature of both dark matter and dark energy are still unknown, yet 
based on their characteristics, different models have been suggested namely tiny positive cosmological constant, 
quintessence, phantoms, chaplygin gas, dark energy in brane- worlds and many others (see the references [4]). 

The dark energy, by virtue of its repulsive gravitational nature, is interesting to study gravitational collapse 

and formation of black hole. As all massive stars do not form black holes (may be neutron stars or white 

dwarfs), so it is generally speculated [5] that the dark energy may play an important role in the collapsing 

H ' stars. 
03 \ 

In the present study, a collapsing spherical star is considered having finite thickness. The star is made of 
dust cloud in the background of dark energy. Let S be the boundary of the star and V + and V~ indicate the 
exterior and interior region of the star. The spherically symmetric inhomogeneous space-time (in V~) can be 
described by the Lemaitre-Tolman-Bondi (LTB) metric element as 

V~ : ds 2 _ = -dt 2 + jr^dr 2 + R 2 dQ 2 (1) 

1 - f(r) 

with R = R(r, t). 

Due to the choice of the co-moving coordinates, the surface £ can be identified as 



E 



a constant. (2) 
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and the metric on it is in the form 

dsl = -dr 2 + RUt) dn 2 (3) 

Thus we have 

t = T and i? s (r) = R(r s ,r) (4) 
In V + , the metric of the space-time in general can be taken as 

V+ : ds\ = -A 2 (T, z)dT 2 + B 2 (T, z)(dz 2 + z 2 dfl 2 ) (5) 
So the bounding surface £ for the coordinates of the exterior space can be expressed as 

z = z (T) (6) 
Then the junction conditions namely dst]^ — cJs+|s give 

■n dT 1 

l) ~dT = i =T= (7) 

ii) R(r s ,t) = Z (T) B(T,z (T)) (8) 

The other junction conditions (relating to extrinsic curvature) due to Israel [6] depend on the choice of the 
space-time outside the star. If the bounding surface S is an energy layer i.e, an infinitely thin matter shell 
appears on S, then the extrinsic curvature is not continuous across S but the jump discontinuity depends on 
the matter of the thin shell over S. On the other hand, the continuity of the extrinsic curvature components 
demand that no such shell exists on S. In fact, whenever the space-time V- is fixed, then the space-time V + 
will determine whether this shell appears or not. In the present paper, it is assumed that no such thin shell 
exists on E and extrinsic curvature is continuous across S. The paper is organized as follows: In Section II, 
basic equations for collapsing inhomogeneous spherical star model are presented. Sections III and IV deal 
with marginally bound and non-marginally bound cases with dark energy in the form of perfect fluid while 
anisotropic fluid form is used for dark energy in Section V. Finally, the paper ends with conclusion in Section VI. 



II. INHOMOGENEOUS SPHERICALLY SYMMETRIC STAR MODEL WITH DARK 
ENERGY IN THE FORM OF PERFECT FLUID 

The energy-momentum tensor for the matter field consisting of an inhomogeneous dust (dark matter) p DM (r, t) 
and homogeneous dark energy in the form of perfect fluid is given by 

Tiiv = (p DM + p + p)u^u v + pg^ (9) 
where in the comoving representation, four- velocity — 5^. 

The explicit form of the energy conservation equations T^. u = for the metric (1) are 

P»M + [^ + R i)p»M=Q(r,t) (10) 



and 
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P+[2^ + ^j(p + p) = -Q(r,t) (11) 

where ( ' ) represents differentiation with respect to time while ( ' ) represents the same with respect to r 
and Q(r, t) stands for the interaction between dark matter and dark energy. 

The non- vanishing components of the Einstein field equations foe the metric (1) having matter field in the 
form (9) are given by (with n = = 1) 

/ f r 2 kk 

V + RR + 1? +2 RR =P ™ +P (12) 
f R? _ R .„ „. 

/' R R 1 R R! . . 

2RR' + ~R~R! + R + R 1 = ~ P ( ' 

In general the interaction term Q(r,t) is non-zero. But if Q(r,t) is zero then from the conservation equation 
(11), R should be in the product form. 

Since in the present work, one mainly considers gravitational collapse, so one assumes R < 0. 
The apparent horizon is characterized by 

R, a R' a = (15) 
If t = t a h(r) be the time of formation of apparent horizon, then from the above equation (15) 

R 2 (r,t ah (r)) = l-/(r) (16) 
The mass function at comoving coordinate V is given by [7] 

m(r, t) = l -R (1 - R, a R' a ) = l -R (r 2 + f(rj) (17) 
So the total mass of the collapsing star at any time r is [8] , 

M(t)= 1 -R^{t)(Rv +/(r)) (18) 

Since the collapsing star is assumed to be not trapped initially (at t = Ti), so on the initial hypersurface 
t = Ti one should have 



R 2 (r,n) + .f(r)-l<0 



(19) 
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III. MARGINALLY BOUND CASE: / = 

In this case the hypersurfaces t ^constant have zero curvature. Now equating the field equations (13) and 
(14), we have a differential equation in R which has a first integral 

_ R R? , . /„„n 
2 r + &=^ (20) 

Integrating once more, the evolution equation for R is 

R , = A*^ rf + 9M (21) 
3 R 

where p(t) and g(r) are arbitrary functions of t and r respectively. 
Now from the field equation (13) using (20), one finds 

p(t) = -fi{t) (22) 

As for dark energy p(t) is negative, so the arbitrary function p(t) is positive. Then from the field equation 
(12) one gets 

p DM +p = p(t) + ^ (23) 

It is to be noted that if a barotropic equation of state of the form p = ep is taken for dark energy, then one 
arrives at a contradiction to satisfy the conservation equations (10) and (11). 

In particular, if p is taken to be constant and the dark energy is in the form of cosmological constant i.c, 
p = —p = n, then from (23) 

P dm - 9'{r)/R 2 R' 

and from conservation equations Q = 0. 

Thus when the dark energy is in the form of cosmological constant then the model corresponds to inhomo- 
geneous dust collapse with a cosmological constant. 

Now, apparent horizon will form if the cubic equation 

pR 3 -3R + 3g = 

has at least one positive root. Various cases are as follows: 

i) For p > and 3g = -^=, a unique apparent horizon will form, given by 

R ah (r) - ^. 

ii) If p > and 3g < -J=, then the cubic equation has two positive real roots, which correspond to 
cosmological and black hole horizons given by 

R c {r) = ^Cos[\Cos-'{-l gs fJi)] 

and 



Mr) - ^Cos [f + \Co8-i (-§Vm)]- 
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iii) For p > and 3g > ^= there are no positive roots and hence there are no apparent horizons. 
The details of this collapsing process has been studied in Ref. [13]. 

The effect on the collapsing star will be studied separately for a dust cloud and dark energy to examine the 
different roles that may play during the evolution. Also their joint effects (with or without interaction) will be 
investigated subsequently. 



Case I: 



(Pdm 7^ 0, p = = p) 



This case corresponds to inhomogeneous dust collapse in LTB model. This has been studied exhaustively by 
Joshi and collaborators [9] (for higher dimensional study, see Ref. [10]). 

Case II: (p DM = 0, p = ep ^ 0) 

As Q(r,t) — [from equation (10)], so area radius R must be in the separable form 



R(r,t)= VW)Rt{t) 
for consistency of the conservation equation (11) and field equations (12)-(14). 
Then from the conservation equation p(t) = p a /R^ 1+€ ^ 
and 



(24) 



p(t) = 



epo 



R 



.3(l+e) 



Thus the differential equation in R t is 



Rt — — R 



e P0 p -(3e+l) , 1 



+ R t 



In the following, this differential equation in R t is solved for e = —1, 0, 1 and — |. 
a) Dark energy in the form of cosmological constant: (e = — 1) 

Equation (26) can be solved easily as 



(25) 



(26) 



= J— Sink 



(t - 1 ) 



or \l — Sink 
Po 



according as R > or < 0. 

The mass function has the expression m(r,t) = \g{r) ^1 + ^p 1 )- 



/3pa 



(to - t) 



(27) 



Hence it represents either an ever expanding model of the universe starting from the big-bang at t — to or a 
collapsing model of the universe having infinite volume at t = — oo and collapse to big crunch at t = to with a 
finite mass. 



b) 



0: 



One can see easily that it corresponds to a pure homogeneous dust (not dark energy) with Rt <~ ta as 
expected in standard FRW model. 



G 



c) Dark energy in the form of perfect fluid: (e = 1) 

Here the scale factor grows with time as 



a? = ? + 7(*-*>) 2 



(28) 



with p(t) = p /R$. 
d)e 

In this case we have 



1 . 

3- 



y/ntfi + Rtpa - Sink- 1 ( v/iZ^) = y (*o - t] (29) 
The mass function has the expression m(r,t) = \g{r) ^1 + £3 ^ ±> j- 
Case III: (p ± 0, p DM ± 0, Q = 0) 

As Q = 0, so similar to Case II, R must be in the product separable form (i.e, eq. (24)) for consistency of 
the conservation equation (10). The expressions for the matter densities become 



Pl 



h(r)Rt 3 and p = p R t 3(1+e) 



(30) 



with po, an arbitrary constant and h(r), an arbitrary function of r alone. Also the evolution equation for R t 
takes the form of equation (26). 

The consistency of the field equation (12) demands p DM to be also homogeneous with h(r) = 3. Thus it 
is very similar to case II and physical consequences are identical. Note that here the role of dark matter is 
insignificant. 

Case IV: (p ? 0, p DM ? 0, Q £ 0) 

For the interacting dark matter and dark energy, the evolution equation for R, the expressions for p(t) and 
(p(t) + p DM ) are given by equations (21)-(23) respectively. 



Now for explicit solution of R one assumes (see Ref. [11]) 

p(t) = /lot", (/io and s are positive constants) 
and the general solution for R has the form [12] 



(31) 



Ri 



' \Tt\cxh 



| S -2| 



2y^ 
1^21 



t ~ 



\s-2\ l 



2V\ 



t 2- 



(32) 



according as £ is an integer, non-integer and s = 2. Hence C\ and C 2 are arbitrary functions of r and 



It can be shown (for details see Ref. [12]) that the model approaches isotropy along the fluid world line as 
t — > oo. 
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IV. NON-MARGINALLY BOUND CASE: / ^ 

In this case, the first integral of equation (13) gives the dynamical equation for the area radius as 



■ 2 p(t)R 2 , 9(r) 



i^ = -i^_ + ^-/(r) (33) 

with g(r) as arbitrary function of r. The explicit form of total energy density is (obtained from equation 
(12)) 

P + p DM = ^j-p(t) (34) 
In particular, if one chooses p(t) = —p(t) , then 

P dm = g'(r)/R 2 R' (35) 

and consequently from the conservation equations (10) and (11) one obtains respectively 

Q = and p ^constant A (say) . 

Hence as in marginally bound case, this can also corresponds to a collapsing dust (inhomogeneous) with a 
cosmological constant. Note that, due to the difference in the evolution equations [see eqns. (33) and (21)], the 
collapsing processes will not be identical. 



One may note that if the separable (product) form of R is assumed, i.c, R = \/ f(r) £(t), then £(i) satisfies 

e + 1 = _ P _^ + vM (36) 

s 

where rj{t) is an arbitrary function of t alone. The form of energy density becomes 

P + P DM = ^-P(t) (37) 

In particular if p(t) = —p(t) = A(i), then 



(DM ^-3 

and the arbitrary function r](t) is related to X(t) by the relation 



(38) 



Thus dark energy behaves as time dependent cosmological constant and the model describes homogeneous 
dust collapse with time dependent cosmological term. 



V. DARK ENERGY IN THE FORM OF ANISOTROPIC FLUID 



Here if the dark energy is in the form of anisotropic fluid with energy momentum tensor 
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T (DE)n = dia 9(-P,Pr,P T ,P T ) (40) 

then the total energy momentum tensor is 

T (I> = Pdm u M u " + T (DE)» ( 41 ) 

So the conservation equation gives 



Pom + (^ + w)pOM=Q(r,t) 



(42) 



and 



R' 

p' r + 2( Pr -p T )—=0 (43) 

One may note that the dark energy in the form of anisotropic fluid must be inhomogcncous, otherwise from 
the conservation equation (43), p' r = implies p r = p T i.e, fluid is isotropic (discussed in the previous sections). 
Conversely, if the dark energy fluid is inhomogeneous (i.e, p' r ^ 0), then again from equation (43), p r ^ p T , i.e, 
the fluid must be anisotropic. 

Now from the field equation (13) (replacing p by p r ), the evolution equation for R becomes (integrating once), 



R -^T~R 



j Pr (r,t)R 2 dR- f(r) (44) 



Now assuming the regularity of the initial radial pressure at the centre and blowing up at the singularity, the 
form of p r can be taken as (for details see Ref. [12]) 

^ ^ ^ (45) 

where a(r) is an arbitrary function of r such that a(r) ~ r" near r = 0. and n(> 0) is any constant. 
Thus the radial velocity of the collapsing shells at a distance r from the centre is given by 

R» = 9 -^-^- n -m (46) 

R 3 — n 

Also from the conservation equations, the expressions for the energy densities become 

_ Po(r) 



+ 



1^ J QR 2 R'dt (47) 



r DM R2RI r2r 

and 

"^r-ikj^ < 48 > 

where H(R,t) — pi(r) — f^-R 3_n and the arbitrary functions po and pi are restricted by the relation 
[consistency of equation (12)] 
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P'l +Pq= g(r) 



(49) 



The expression for tangential stress becomes 



a(r) 



1 



nE! 



+ 



(n ± 3) 



(50) 



Pt = 



R n 



2(R' + Rv') 



2R n ~ 1 {R' + Rv') ' 



A similar collapsing process has been studied extensively in Ref. [12] for marginally bound case only (/ = 0) 
and it is calculated that in general, pressure tries to resist the formation of naked singularity. Thus formation 
of a black hole from the collapsing star due to the presence of dark energy in the form of anisotropic fluid is 
more favorable than formation of a naked singularity. 



The paper deals with a detail study of an inhomogeneous spherically symmetric star model having dark 
matter in the background of dark energy as the matter content. The dark matter is taken in the form of 
inhomogeneous dust while for dark energy, both homogeneous perfect fluid and anisotropic fluid models are 
considered separately. One may note that the present study can easily be extended to quasi-spherical collapsing 
star models and the results will be identical. 

For dark energy in the form of perfect fluid, both marginally bound (/ = 0) and non-marginally bound 
(/ 7^ 0) cases are considered in two different sections. When / = 0, dark matter and dark energy are considered 
both separately and in combined form. It is observed that due to the presence of dark energy, trapped surfaces 
do not form at all and collapsing models lead to big crunch singularity. In a particular situation (when 
e = —1), the non-marginally bound case represents inhomogeneous dust collapse with a cosmological constant 
while assuming separable form of area radius, the model corresponds to homogeneous dust collapse with time 
dependent cosmological term. As for inhomogeneous dust with cosmological constant, trapped surfaces may be 
possible while for homogeneous dust with time dependent cosmological term, formation of trapped surfaces is 
not possible. Hence naked singularity is more favorable in second case than in first case. 

Lastly, when anisotropic fluid represents dark energy, it is found that black hole formation is more favorable 
than naked singularity. Therefore, one may conclude that inhomogeneity, both in dark matter and dark energy 
favors formation of trapped surfaces, while homogeneous cases support the singularity to be naked. 
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